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Bilayer graphene as an helical quantum Hall ferromagnet 
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The two-dimensional electron gas in a bilayer graphene in the Bernal stacking supports a variety 
of uniform broken-symmetry ground states in Landau level A'' = at integer filling factors v £ 
[—3, 4] . When an electric potential difference (or bias) is applied between the layers at filling factors 
1/ — —1,3, the ground state evolves from an interlayer coherent state at small bias to a state with 
orbital coherence at higher bias where electric dipoles associated with the orbital pseudospins order 
spontaneously in the plane of the layers. In this paper, we show that by further increasing the 
bicis at these two filling factors, the two-dimensional electron gas goes first through a Skyrmion 
crystal state and then into an helical state where the pseudospins rotate in space. The pseudospin 
textures in both the Skyrmion and helical states are due to the presence of a Dzyaloshinskii-Moriya 
interaction in the effective pseudospin Hamiltonian when orbital coherence is present in the ground 
state. We study in detail the electronic structure of the helical and Skyrmion crystal states as well 
as their collective excitations and then compute their electromagnetic absorption. 

PACS numbers: 73.21.-b,73.22.Gk,72.80.Vp 
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I. INTRODUCTION 

Experiments^! have shown that, when placed in a per- 
pendicular magnetic field, a graphene bilayer supports 
a set of Landau levels with energies given by E^ = 
±nuJc\/\N\{\N\ + l) where iV = 0, ±1, ±2, . . . All lev- 
els except iV = are fourfold degenerate in addition to 
the usual degeneracy related to the guiding center coordi- 
nate. An electronic state is thus specified by its Landau 
level index N, its guiding center index X (in the Lan- 
dau gauge), its spin index a — ±1 and its valley index 
^ — zLK. In N = 0, we must add an additional orbital 
index because states with both n = and n = 1 Landau- 
level character have zero kinetic energy. The eightfold 
degeneracy of the N = Landau level is revealed ex- 
perimentally by a jump in the quantized Hall conductiv- 
ity from —4 (e^/ft.) to +4 (e^/ft.) when the charge den- 
sity is tuned across neutrality in moderately disordered 
samples?. 

Some of usS have shown that the close proximity be- 
tween the two graphene layers in bilayer graphene leads 
to spontaneous interlayer coherence in iV = when 
Coulomb interaction is taken into account. In a pseu- 
dospin language where the spin, layer and orbital degrees 
of freedom are each mapped to a 5^ = 1/2 pseudospin 
{i = tr, ^,n), electron-electron interactions at integer fill- 
ing factors completely lifts the degeneracy of the bilayer 
octet producing spontaneously broken-symmetry ground 
states with spin, valley and orbital polarizations. In con- 
sequence, quantum Hall plateaus should occur at all in- 
teger values of the filling factors from v = —3 to ly = 4 
in Landau level iV = 0. The existence of these additional 
plateaus has recently been confirmed experimentally in 
suspended bilayer graphene samples and bilayer graphene 
on Si02/Si substrates^. 

The possibility to study novel broken-symmetry (BS) 
states is always exciting and the new BS states in bi- 
layer graphene are no exceptions. For example, the 



interlayer-coherent state at z/ = — 3 (and i^ = 1 if the two- 
dimensional electron gas (2DEG) can be considered as 
fully spin polarized) is a quantum Hall layer-pseudospin 
ferromagnet with pseudospin wave dispersion lo ^ q^ 
that contrasts with the usual linear dispersion found 
in interlayer-coherent states in semiconductor double- 
quantum- well systems. This unusual dispersion can be 
related^ to a vanishing of the counterflow superfuid den- 
sity. At filling factor v — —1,3, the interlayer-coherent 
state has the usual linear dispersion^. If a positive elec- 
tric potential difference A^ (which we refer to, in this 
paper, as the bias potential) is applied between the lay- 
ers, the charge is progressively transferred into the bot- 
tom layer. At i/ = —3, 1 the ground state with all the 
charge in one layer supports an orbital pseudospin mode 
that can be viewed as an intralayer cyclotron resonance^. 
This mode is gapped due to the finite bias and should be 
detectable in microwave absorption experiments. 

The broken-symmetry states related to the orbital de- 
gree of freedom are especially interesting. For a spin- 
polarized 2DEG, they occur at filling factors zy = —1,3 
and above a certain critical bias A^ where one of the 
layer is completely filled. The homogeneous orbital 
broken-symmetry states have a finite density of electric 
dipoles that collectively order in the x — y plane^i^. The 
orbital pseudospin mode corresponding to the collective 
motion of these dipoles is gapless despite the finite bias. 
It is a Goldstone mode due to the breaking of the U{1) 
symmetry of the pseudospins in the x — y plane. 

It was shown recently^, that there exists a 
Dzyaloshinskii-Moriya (DM) interaction between 
the orbital pseudospins at v — —1,3 that causes the 
orbital pseudospin mode to soften at a finite wave 
vector q as the bias A^ is increased. It was conjectured 
that, above a critical bias, the ground state should 
be some kind of helical state. In the present work, 
we find a more complex scenario. Working in the 
Hartree-Fock approximation, we find that, as the bias 



is increased, the ground state follows the sequence of 
transitions: uniform state - Skyrmion crystal - helical 
state - Skyrmion crystal - uniform state. The phase 
diagram is symmetrical about the bias A^ /2 where 
the charge is equally distributed in both orbitals n — 
and n — 1. Interestingly, our phase diagram looks very 
similar to that found recently in a thin film of the helical 
magnet Feo.5Coo.5Si where a Skyrmion crystal as well 
as an helical and a ferromagnetic uniform phases have 
all been observed using Lorentz transmission electron 
microscopyZiS.. In this system, the phase transitions are 
induced by a transverse magnetic field. In the bilayer, 
the role of the magnetic field is played by the bias. The 
effective pseudospin hamiltonian of the orbital states in 
a graphene bilayer is quite complex. Since charge and 
pseudospin are coupled, a pseudospin texture such as 
that found in the orbital Skyrmion crystal produces a 
charge density in real space and so the direct (or Hartree) 
Coulomb interaction must be considered together with 
the other competing interactions mentionned above. 

In this work, we study in detail the helical and 
Skyrmion crystal states. We derive their electronic prop- 
erties as well as their collective excitations and compute 
their electromagnetic absorption. We show that the effec- 
tive pseudospin model that describes these states involves 
nonlocal as well as anisotropic exchange interactions be- 
tween the pseudospins. These Coulomb exchange inter- 
action tend to align the orbital pseudospins while the 
DM term favors a rotation of the pseudospins in space. 
The resulting ground states result from a balance be- 
tween these competing forces just as in helical magnets 
such as MnSi and Fei_j,Coj:Si. One major difference in 
the graphene bilayer is that the DM interaction in the 
orbitally coherent state is of an entirely different origin 
than the spin-orbit interaction at work in usual helical 
magnetsi". 

Our paper is organized in the following way. In Sec. 
II, we study the non-interacting states of the graphene 
bilayer within a two-band low-energy model. In Sec. Ill, 
we discuss the validity of our approximation of assum- 
ing complete spin polarization. This discussion is needed 
since increasing the bias pushes down (up) half of the spin 
down(up) levels. Some levels of opposite spin cross at fi- 
nite bias and this can introduce new phases in the phase 
diagram or make some phases disappear—. In Sec. IV, we 
summarize the Hartree- Fock approximation as well as the 
generalized random-phase approximation (GRPA) which 
we use to study the collective excitations. In Sec. V, we 
study the electronic properties of the different phases as 
well as their collective excitations. We discuss their elec- 
tromagnetic absorption in Sec. VI and conclude in Sec. 
VII. 



II. EFFECTIVE TWO-BAND MODEL 

We consider a graphene bilayer in Bernal stacking as 
represented in Fig. [TJ The bilayer is placed in an ex- 



ternal transverse electric field in order to allow an elec- 
trical potential difference As between the layers and in 
a transverse magnetic field B. We denote the two basis 
atoms in the top layer by Ai and Bi and those of the 
bottom layer by A2 and B2 with atoms Ai sitting di- 
rectly above atoms B2 as shown in the figure. The band 
structure of this system is calculated in the tight-binding 
approximation with the hopping parameters: 70 = 3.12 
eV,7i = 0.39 eV,74 = 0.12 eV,A = 0.0156 eV, taken 
from Ref. lid . The parameter 70 is the intralayer hopping 
term between nearest-neighbor carbon atoms, 71 is the 
interlayer hopping between carbon atoms that are part 
of a dimer, 74 is an interlayer hopping term between two 
carbons atoms that are not part of a dimer {Ai — A2 and 
Bi — B2) and A represents the difference in the crystal 
field experienced by the inequivalent atoms Ai and Bi in 
the same plane. We neglect the trigonal warping term 
73 (the A2 — Bi hopping), a correct approximation at 
sufficiently high magnetic fieldii. 




FIG. 1; (Color online) Lattice structure of the bilayer 
graphene in the Bernal stacking. The spheres represent car- 
bon atoms. 

The electronic dispersion consists of four bands. In 
the absence of bias and magnetic field, two of these 
bands meet at the six valley points K,K' of the recip- 
rocal lattice. Here, we take as the two non-equivalent 
points K = (-47r/3ao,0) and K' = -K = (47r/3ao,0) 
where ao = 2.46 A is the lattice parameter of graphene. 
The dispersion of the two low-energy bands, for small 
wavevector k measured from either K or K', is given by 
E^v^ (k) = ±fi2/c2/2m* with valley index ^ = ±1. The ef- 
fective electronic mass is defined by m* = 2fi'^ji/3jQaQ = 
0.054mo where itiq is the bare electronic mass. The two 
high-energy bands are separated from the two low-energy 
bands by a gap of order 71 . 

To describe the low-energy {E << 71 ) excitations of 
the tight-binding model, we assume complete spin polar- 



ization and use an effective two-band model^'' where the 
Hamihonian is given by 



H^K 



-^^ + (C + e/3AB)aat 



e^ + (C-e/3As)a^a 

(1) 
in the basis (^2,-81) for H^ and (Si, A2) for -ff°K- In 

Eq. ([T]), a,a'^ are the ladder operators and we have de- 
fined the parameters 



C = 2sgn(7o74)V/3/347i+/3A, 

where sgn denotes the signum function and 

fjfj* 
/3 = — ^ =8.86x IQ-^B, 

71 

p^ = f^iVM^ 1.31x10-^5, 
V70/ 71 



(2) 

(3) 
(4) 



(with B in Tesla). The effective cyclotron frequency is 
defined by uj* = eB/m*c. 

When 74 = A = A^ = 0, the Landau level energies in 
each valley are given by 



Ei^=±./\N\m 



i)f^: 



(5) 



with N — 0,±1,±2, ... where N is the Landau level. If 
the Zeenian coupling is taken as zero, all Landau levels 
are thus four time degenerate (including spin and valley 
degrees of freedom) with the exception of iV = that 
is eight times degenerate as shown below. With finite 
values of 74, A or A^, the degeneracy is lifted and we 
find for the states in A'^ = the following spinors and 
energies (we use the Landau gauge with A = (0, Bx, 0)): 



vZ(r))'^«Vo.x,. - l^As-aA^ 



(6) 



( h,.Z (r) ) ' ^"-'^^-- = b""^ - ^^^ - ^^^^ + ^^ 

where we have added a Zeeman coupling A^ = 
gHBB/2 = 0.58 x lO^^S eV (with B in Tesla) for more 
generality. In the absence of couplings and with ^ = 0, 
the A^ = Landau level has an extra orbital degeneracy 
since the two spinors above are then degenerate. The 
wave functions hn,x (r) = e~^^yl^ ip^ {x — X) /a/L^ are 
the eigenstates of an electron with guiding center X in 
the Landau gauge, and ipn (x) is the wave function of 
the one-dimensional harmonic oscillator. Note that with 
our choice of normalization for the functions Lp^ (x) , the 
action of the ladder operators on the states tpn (x) is 
given by a^n (x) = W^ + I'/'n-i-i (x) and aipn (x) — 

-iy/nipn-l {x) . 

We only consider Landau level A'^ = in our work 
so that, from now on, we will drop the index N. The 
index n will always refer to the orbital states n = 0, 1 
i.e. to (po (x) and ipi (x) . Note that it is clear from Eq. 
([5]) that the valley K(K') eigenstates are localized in the 
top(bottom) layer. For iV = 0, the layer index is thus 
equivalent to the valley index. 



III. QUANTUM HALL FERROMAGNETS 

To describe the electronic phases in Landau level N — 

we use a pseudospin language where we associate to the 
layer ^ = ±K (or valley) and orbitals n = 0, 1 a spin-half 
pseudospin. When Coulomb interaction is considered in 
the Hartree-Fock approximation, the ground states at 
zero bias and C = obey a set of Hund's rules in which 
spin polarization is maximized first, then layer polariza- 
tion to the greatest extent possible, and finally orbital 
polarization to the extent allowed by the first two rules^i. 
In the pseudospin language, the ground states are thus 
various types of quantum Hall ferromagnets (QHF's). 

At zero bias and with C = 0, interlayer coherence is 
present in the ground state at all integer filling factors 
but there is no orbital or spin coherence. Interlayer co- 
herence occurs naturally because of the proximity of the 
two layers (the interlayer separation d = 3.34 A is very 
small in bilayer graphene) but disappears quickly when 
a finite bias A^ is applied. Adding a finite bias A^ en- 
riches considerably the phase diagram with new states 
such as orbital and spin QHF's. 

In previous works^, it was found that above a criti- 
cal bias Ab {i'), the ground state is an orbital QHF with 
Ising character (i.e. no orbital coherence) for ly — —3 
and a ground state with an a; — y character (or orbital 
coherence) for j/ = — 1 . The dispersion relation of the or- 
bital pseudospin mode is gapped at w (q = 0) = /3Ab in 
the former case and gapless with an highly anisotropic 
dispersion at small wave vector in the latter case. To 
get these results, it was assumed that complete spin po- 
larization holds even at finite bias and the correction ^ 
was not considered in the single-particle energies. When 
these two assumptions are relaxed, the phase diagram at 
Ab 7^ is modified. An exhaustive study of the phase 
diagram for uniform states at all integer filling factors 
ly E [^3,4] will be presented elsewhere^. We want to 
mention here, however, some changes that occur in the 
phase diagram. 

First of all, we remark that orbital coherence is driven 
by the /3Ab term in Eq. ^. Together with the term 
C, they lift the degeneracy between levels n = and 
n — l.The ordering of the non-iteracting levels is then as 
illustrated in Fig. [H 

When the correction C, is neglected, the orbital splitting 
/3Ab << Ab- For positive bias, the n — orbital state 
in the bottom(top) layer is lower(higher) in energy than 
the n = 1 orbital state. The orbital energy splitting has 
thus different signs in the K and K' valleys. No orbital 
coherence is possible at i/ = — 3 because level n = 1 is 
always above level n = at all bias. In that case, the 
ground state has all electrons in valley K' with n = 
because the Coulomb exchange energy is more negative 
in 71 = than in n = 1. The situation is different for 
1/ = — 1. In that case, the valley K' is filled above a 
critical bias and the remaining electrons occupy valley K 
where level n = 1 is now below n = 0. Now, because of 
the better exchange interaction in n = 0, the eigenstates 
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plete spin polarization. Our results should apply with 
some minor changes to the orbital states at z^ = 1 and 

In concluding this section, we remark that the phase 
diagram for filling factors v — —1 has been studied 
in some details before^;^ but only the homogeneous 
states have been considered. Non-uniform states such as 
Skyrmion crystals with valley or orbital pseudospin tex- 
tures have also been considered but near integer filling 
factors onlyi^. In this work, our focus is on the non- 
uniform states at precisely v = —I. 



IV. FORMALISM 

In the rest of this paper, we concentrate on the study 
of the orbital coherent states at z^ = — 1 , assuming C = 
and complete spin polarization. Our calculations are all 
done at zero temperature. 



FIG. 2: Energy levels for one spin state in the presence of a 
finite bias. 



in K are bonding and anti-bonding states of n = and 
n — 1 with the electrons mostly in n = at low bias and 
mostly in n = 1 at large bias. This produces an orbital 
coherent state. 

If C 7^ 0, we see from Eq. © that the only effect is 
to push the orbital coherent states to higher bias. In 
itself, this is not dramatic provided our effective two- 
band model is still valid at the new critical bias. But, if 
the bias is increased, the spin degree of freedom must 
be considered. At zero bias, the four spin up states 
are separated in energy from the four spin down states 
by an exchange-enhanced Zeeman gap of order (? j kI 
where k is the dielectric constant of the substrate and 
i = yJhEJeB is the magnetic length. With bias, two lev- 
els with spin up(down) are shifted upward(downward) 
in energy. When levels with different spins cross, states 
with spin coherence [x — y spin QHF) become possible 
and they may replace the orbital state in the phase dia- 
gram. Whether this is the case or not must be decided 
by numerical calculations. For 74 — 0.12, a numerical 
calculation^ shows that the orbital coherent state is ab- 
sent for V = —\. We remark, however, that the exact 
value of 74 is not precisely known and than orbital co- 
herence does occur at smaller values of this parameter. 

Fortunately, the phase diagram for the eight-level sys- 
tem is very rich and numerical calculations^ show that 
orbital coherence with no interlayer or spin coherences 
occurs in some range of bias but at the higher filling 
factors V = \ and v — "i and even when 74 is as large 
as 74 = 0.12 and spin mixing is considered. Since we 
are confident that orbital coherent states do occur in the 
phase diagram in bilayer graphene, we will study, in this 
paper, the simplest case with z/ = — 1,74 = and com- 



A. Hartree-Fock Hamiltonian 

With the restrictions outlined above, the analysis is 
reduced to that of a two-level system in valley K since 
the two filled levels in valley K' can be considered as 
inert. The non-interacting Hamiltonian is given by Eq. 
([1]) while the many-body Hamiltonian in the Hartree- 
Fock approximation is^ 

HhF = N^EnPn^r (0) (7) 

+^vX!-^"i'"2,n3,"4 (q) {Pni.n2 (-Q)) Pn3,«4 (q) 



-iV^^X„,, 



(q)(pni,«2 (-q)> p«3,n4 (q) : 



where n^ = 0, 1 is the orbital index and N^p = S/2tt£'^ is 
the Landau level degeneracy. The single-particle energies 
can be simplified to En = —nfSAs- In deriving Eq. (O, 
we have taken into account a neutralizing positive back- 
ground so that the q = contribution is absent in the 
Hartree term. This convention is indicated by the bar 
over the summation. 

The density operators in Eq. ([7]), are defined by 



'"^' (q) - M E_ 



e 2 



Iq^Jfi+Xs) 



(8) 



Xi,X2 



t s: 

'CXi,niCX2,n2'5Xi,X2+g„f2, 



where c\^ „ (cx.n) creates(destroys) an electron in state 
{X, n) in the Landau gauge. The Hartree and Fock in- 
teractions are given by 



iif«i,«2,n3,«4 (q) = ( ^ ) -^-ft^ni,n2(q)-ft^«3,n4(-q)(9) 

X (a) - f ^^^ H ('D)e''>''P^'ClO^ 

^^ni,n2,n3,"4 Vm/ / -'-'"l ,"2,n3,n4 vfy '^ X'-^J 



where k is the effective dielectric constant at the position 
of the graphene layers. The Coulomb energy e'^ /k£ = 1 
1.3a/B mcV with B in Tesla and we take n ~ 5. 

The form factors which appear in H and X are given 
by 



^0,0 (q) = exp 



2«2 



-q'e 



(11) 



Ki,i (q) = exp I — ■ — I ( 1 



2 ; ' 

^1,0 (q) = I ^ ) exp 



V2 
Kqs (q) = I —7^ ) exp 



V2 



They capture the character of the two different orbital 
states and play an important role in the physics of the 
orbital phase. Detailed expressions for the Hartree and 
Fock interactions H and X are given in Appendix A of 
Ref. i. 



B. Order parameters in the coherent phases 

The order parameters of the orbital phases are ob- 
tained from the single-particle Green's function 

G„,„. {X, X', t) = - (TrCn,x (r) cl,,, (0)) , (12) 

where Tt is the imaginary time ordering operator. We 
define the Fourier transform of the single-particle Green's 
function as 



G„x(q,T) = ^ E ^"^"^''^'''^ 



(13) 



X,X' 



^5x.X'-qypGn,n' {X,X',t) 



SO that the order parameters of the coherent phases are 
related to the Green's function by 



{Pn,n' (q)) = Gn',n (q,T = 0" 



(14) 



The equation of motion for the Green's function in 
the Matsubara formalism and in the Hartree-Fock 
approximation is given by 



(iw„ -I- l-l/h) Gn,n' (q, i^n) " (15) 

E E "^"'™ (^^' ^') 7q,q'G'm,„' (q, iw„) == (5„,„'(5q,0, 
m— 0,1 q' 

with the matrix elements 

T„,™(q,q') = <„(q-q')-<™(q-q') (16) 

— plS.B5c^,c^'5n,m5n.l, 



and the phase factor 



7q,. 



-jqXq'£^/2 



(17) 



The Hartree and Fock potentials are defined by 

Un:,.nA'^) = E -^"i'"2,"3 ,"4 (-q) (pm ,"2 (q)) ,(18) 

ni,"2 
C^n3,n4(q) = E ^"I'»4."3,n2(-q)(pni,n2(q))-(19) 



The self-consistent Eq. ([15]) must be solved numeri- 
cally in an iterative way in order to get the order param- 
eters in the different orbital phases. Once this is done, 
the Hartree-Fock energy is obtained from 



E 



HF 



N 



-/3Ab (pi,i (0)) 



^E E ^«. 

q ni ,...,n4 



(20) 

(q)(Pni,n2(-q))(Pn3,n4(q)) 



~ 2 ^ ^ ^«i,n4,«3,n2 (q) {Pnun2 (-q)) {Pns-^u (q)) • 
q ni,...,n4 

The Hartree-Fock equation of motion for the Green's 
function leads to the sum rule (at T = OK): 

EEK^"-"^(q))i' = (^"-rn(o)), (21) 

q "2 
and we have by definition 

{pn,n (0)) = Vn, (22) 

where Vn is the filling factor of the n level. 

C. Induced electric dipoles 

In the pseudospin language, we associate state n = 
with pseudospin up and n = 1 with pseudospin down. 
We then define the guiding-center density operator 



p«(q) = po,o(q) + pi,i(q) 



(23) 



and the three components of the pseudospin density vec- 
tor p (q) are given by 

Pxiq.) = ^{po,ii<l) + Pifi{cL)), (24) 

Py(q) = ^ (po,i (q) - pi,o (q)) , 
Pzin) = 2 (^0,0 (q) - pi,i (q)) • 

The coupling of the 2DEG with a uniform external 
electric field in the plane of the layers is given by He = 
— e / dm (r) (r) with E — —Wcj) (r) . In this expression, 



we must use the true density n (r) which is the Fourier 
transform of 

n (q) ^ N^Y^ p„,.n2 (q) Kn,,n2 (-q) • (25) 

ni,n2 

This coupling can be written as 



H = -eN^J dvp,Ar)cj,{r) (26) 

+V2£eN^ {p., (q ^ 0) E^ ~ py (q = 0) Ey) , 

where we have defined p^ (q) — exp (— (7^^^/4) pi (q) with 
i = n,x,y,z. In the three states studied in this pa- 
per (uniform, hehcal, and Skyrme crystal), the guiding- 
center density (p„ (r)) is uniform and so (p„ (q)) = Sq.Q- 
For these states, the couphng with the electric field is 
simply given by 

He = -N^d (0) • E, (27) 

where we have defined the electric dipole operator 

d (q) = -~eV2ee^^ {p, (q) , -py (q)) . (28) 

The fact that orbital coherence leads to a finite density of 
electric dipoles in the plane of the layers was first shown 
in Ref. M- 

Note that for the states with (p„ (q)) ~ 6qfi, the sum 
rules of Eq. ([2T|) are, in pseudospin language, equivalent 
to 



Xni,n2,n3,n4 (qj q i *"nj 
Xni,n2,n3,n4 vQi *1 ' *"nj 



(31) 



n 2^ /^Xrii,n2,"6,n6 ('^''^ ;"^n) 



IB, •••,"8 q" 



xH, 



(q")x« 



(q",qVaO 



"fi XI H^"i."2,«5,n6(q'q";*^") 
"5, •••,"8 q" 

X^n5,T!8,n7,"6 vQ / Xn7,"8,n3,"4 Iq : q I *^^rij j 

where ri„ is a bosonic Matsubura frequency and 
the Hartree-Fock two-particle Green's function 
X°i,«2,n3,«4 (q>q';*^n) is given by 



[inn,, - {Er,, - E,,,)]xt,n2,n3.nA^^^' ^^n) (32) 

?l(/3ni,n4(q-q'))^n2,«37q,q' 

-ft (Pn3,n2 (q " q')) ^ni,n4 7q,q' 

-T.T.u"n, (q- q")7q,q"X°,„2,n3,n4 (q">q'>^n) 

n q" 

+ Yl Yj^Ln (q - q") 7q,q"x°i,n,„3,n4 (q"' q'' ^") 

n q" 

+ E E <ni (q - q") 7q.q"X°.„2,«3.n4 (q"' q' ' "") 
n q" 

- E E ^'^2,n (q - q") 7q,q"x"i,„,„3,n4 (q"' q'' ^") ■ 



EK'^-(q))i' + i('^«(q))i' + i('°^(q))i' = i 



(29) 



which is an average normalization condition for the pseu- 
dospin field vector i.e. / dv {p (r)) = 1/45' where S is 
the area of the 2DEG. The modulus of the pseudospin 
field is, in general, not constant in space. 



Note that the response functions depend only on the or- 
der parameters {pn,m (q)) computed in the HFA. Eqs. 
pil32p can be solved numerically by writing them in 
a matrix form. The collective excitations are then 
given by the poles of the retarded Green's function 
Xni,n2,n3,n4 (q'q''"^) ^hich is obtained by the analytic 
continuation ii7„ — > uj + iS oi the corresponding two- 
particle Green's function. 



D. Collective modes 

To study the collective excitations, we compute the 
two-particle Green's function 



x«i,"2,n3,"4 (qiq';'r) (30) 

= -N^ {TrPm.n2 (q,T)/o«3,«4 (-q'.o)) 

+ N^ (Pm.na (q)) (P«3,ru ("q')) 

in the generalized random-phase approximation (GRPA). 
In this approximation, X"i,n2.n3,"4 (q^q's''') is the solu- 
tion of the equation 



V. PHASE DIAGRAM FOR THE ORBITAL 
PHASES AT 1/ = -1 

We define an average pseudospin field by 

p (q) = ^/2fe (- {p., (q)) , {py (q)) , (p, (q))) , (33) 

so that the dipole field can be related to the parallel 
component of this vector by 



d (q) = e '4 p|| (q) 



(34) 



After a lengthy calculation, we can put the Hartree-Fock 
energy per electron given by Eq. (|20p in the pseudospin 



form 



Ehf ,(e^ 

N ' \Ki 



11 ^ 1 
+i/? (ab - ^A^'^'l p^ (0) 



(35) 



2a 



i^EpiK-q)- [«('?)! +^('?)^(q)]-pii(q) 



7r^^c{q)p^{-q)p,^{q) 



2a 

i 
'4^ 



J2 d {q) (2 X q) • (p (-q) X p (q)) , 



In this equation, the bias As — As/ (e^/k (£)) , a ^ 
\/2ie, I is the 2x2 unit matrix and A (q) is given by 



A(q) 



_ f cos {2ipq) sin {2ipq) 

sin (2(/5q) — cos (2(/3q) 



(36) 



where (fq is the angle between the wave vector q and the 

X axis. The bias A^ /2 is defined as the bias for which 
the charge is equally distributed between the two levels 
n = 0, 1 in the uniform phase. Finally, the pseudospin 
interactions are given by 



a{q) = 2 (iJo, 14,0 (g) - ^1,1:0,0 (g)) 



(37) 



q£e 
-2 






b{q) = 2e2''^- (iJi,o4,o (q) - ^1,0,1,0 (q)) (38) 

dyy^e-y"/^J2 [qty) , 



c{q) 



^ q£e 2 - 



-ffo. 0,0,0 (9) — -'^o, 0,0,0 (q) 

+i?i, 1,1, 1(9) -^1,1,1, 1(9) 

-2 (i/1,1,0.0 (9) - Xo,i.i,o (9)) 
q3£3 ,2,2 



(39) 



d{q) 



J dy(2-2y' + ^je-y'/^Jo{qey), 

-Aie'^^ (-ffi,o,o,o (q) - ^1,0,0,0 (q)) 

-4ze"^- (ffi,i,i,o (q) - ^1,1,1,0 (-q)) (40) 



4 q^e _,2. 
-e 2 



V2 2 




FIG. 3: (Color online) Effective pseudospin interactions as a 
function of wave vector qi. 



These interactions are plotted in Fig. [31 Their values at 
q = 0area(0) = -/f,6(0)==0,c(0)==-fv/f,d(0) = 
0. 

Eq. (PSj) is the energy functional of an effective helical 
pseudospin ferromagnet. Apart from the constant terms, 
there are four distinct contributions to the total energy: 

1. The term /3 TAs - A^ V2) P^ (0) is an effective 

Zeeman coupling that changes sign at the bias 

Ag 72. For As < A^ /2, the pseudospins order 
along the positive z axis. They order along ~z for 

AB>Ai'V2. 

2. The terms involving a{q) and c{q) are non lo- 
cal Hcisenberg exchange interactions between the 
pseudospins. 

3. The term with h (q) A (q) is a dipolar interaction. 
In the small wave vector limit, only the Hartree 

,2f2 

term q£e 2— in Eq. ((38)) contributes significa- 
tively to b {q). To first order in q£, b (g) — >■ qi. The 
dipolar term can be related to the Fourier trans- 
form of the dipole-dipole electrostatic interaction. 

4. The fourth term is a Dzyaloshinskii-Moriya (DM) 
interaction between the pseudospins. With the re- 
definition py — )■ —px and px -^ Py, it takes the 
more conventional DM form which is the Fourier 
transform of I? J dr (p • (V x p)) where D is a con- 
stant. In our case, the DM constant D becomes a 
non local function d (r — r') . It is interesting to re- 
mark that a DM occurs in our model although we 



are not dealing with real spins or spin-orbit interac- 
tion. The physical origin of this term was discussed 
in Ref. |5l From a microscopic point of view, d (q) 
contains only the interactions i?ni,«2,n3,ra4 (q) and 
^rai,n2, "3,714 (l) that do not conserve the orbital in- 
dex (for example: Xi^o.o.o (q))- These terms arises 
in our system because we are dealing with two dif- 
ferent orbitals ho (r) and hi (r) . On the contrary, 
the other effective interactions a{q) ,b (q) and c (q) 
consist of terms that conserve the orbital index. 

The exchange interactions tends to align the pseu- 
dospins together while the DM term favors a rotation 
of the pseudospins in real space. This type of competi- 
tion is usually present in helical magnets such as MnSi 
and Fei_a;Coa;Si. In fact, our numerical calculation gives 
a phase diagram which is similar to that found recently in 
the helical magnet Feo.5Coo.5Si where Skyrmion crystal 
and an helical phase have been observed using Lorentz 
transmission electron microscopy?. 

Because the instability of the uniform phase occurs at 
qy£ w —2 (see below), it is not useful to derive a long- 
wavelength approximation of Eq. (|35|) . The calculation 
of the optimal wave vector q£ for the helical state must 
be done numerically. To obtain the phase diagram for 
the orbital phases, we evaluate the Hartree-Fock energy 
for (once again, (p„ (q)) = Sq^o in all three cases): 

1. The uniform phase (UP) defined by p (q) — po^^q.o 

2. The Skyrmion crystal phase (SKP) defined by the 
set of order parameters {p (G)} where G is a recip- 
rocal lattice vector of the crystal. The triangular 
crystal has lower energy than the square or rect- 
angular lattices. (We have not tried other lattice 
types however.) 

3. The helical phase (HP) defined the order parame- 
ters {p (q)} where q = nqo where n = 0, ±1, ±2, ... 

Within this set of states, we find the following ordering: 



As < 0.06 


UP 


0.06 < As < 0.52 


SKP 


0.52 < As < 3.02 


HP 


2.98 < As < 3.44 


SKP 


3.44 < As < 3.536 


UP 



Above As = A^ = 3.536, all electrons are in level 
n = 1 and there is no orbital coherence. We remark 
that the phase diagram is symmetrical with respect to 

the bias A^ /2 (see Eq. (|43|) below). When the layer 
index is added to the picture, the ground state is an 
interlayer coherent phase with no orbital coherence for 
Ab < 0.0022. Also, although the guiding-center den- 
sity (p (r)) is uniform in space, the real density n (r) is 
not. For (p (r)) to be constant in space, we must have 
(Po,o (q)) = — (Pi,i (q)) for q 7^ in the Skyrmion and 
helical phases. 



A. Uniform phase 

The uniform phase has all orbital pseudospins pointing 
in the same direction in space. Its energy is given by 



Ehf ,(e^ 



(41) 



11 
32 
1 
a 
1 
'2^ 



1 



/3 A 



1^(1) 



P.(0) 



2 P||(0)-p||(0) + -p,(0)p,(0) 



The Hartree-Fock solution for the order parameters is 

(42) 



Pz{^) = 9 -=71)' 



|P||(0)| 
At the critical field 



J-(P.(0))'. 






(43) 



all electrons are pushed in level n = 1 and the orbital 



coherence is lost. For B — IQ T, we find that A 
3.536. From Eq. (|42l) . we have the symmetry 



■(1) 



|P||(0)|a. 
Pz (0)a, 



|P||(0)| 

'Pz (0) 



A^„^'-A, 



A^„^'-A, 



(44) 

(45) 



The band structure consists of two dispersionless bands 
with energy 



E+ 



TT e 
2 V^ 



The energy in the middle of these two bands is 



1 



(E, 



E) 



We obviously have the symmetry 



E. 



2A 



(*) 



£;_. 



(46) 

(47) 

(48) 
(49) 



As can be seen from Eq. (pij) . the energy of the UP 
is independent of the orientation of p|| (0) in the x — y 
plane. The UP is a quantum Hall orbital pseudomagnet. 
It follows that this phase supports a gapless x — y orbital 
pseudospin wave excitation, a Goldstone mode related to 
the breaking of the rotation symmetry about the z axis. 



The pseudospins remain parallel during their motion in 
this mode (for q = 0) and so the DM term does not 
open a gap as it does for the corresponding mode in the 
helical state (see below). The contribution of the DM 
is strongest at larger bias and shorter wavelength. The 
mode's dispersion relation is plotted in Fig. |4] for dif- 
ferent values of the bias. It is highly anisotropic^ with 
an unusual square root dispersion in the qy direction if 
the dipoles are oriented in the direction — x. The DM 
instability of the UP occurs at A^ = 0.37 (as well as 

at Ag' — 0.37 if the bias is decreased from A^') and at 
finite wave vector qyi w —2. This suggests a transition 
to a ground state where the orbital pseudospin field is no 
longer uniform as in the HP. In reality, we find that this 
instability is preempted by a transition to a Skyrmion 
crystal at A^ = 0.06 as we discussed in the previous 
section. 
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FIG. 4: (Color online) Dispersion relation of the orbital pseu- 
dospin mode in the uniform phase for different values of the 
bias. For (1) to (3) respectively: As/ {e^/ntj = 0.1,0.3,0.37. 



where ei,e2 are orthogonal unit vectors and a is a nor- 
malization factor. The last term with rj describes the 
uniform polarization due to the bias. It makes the mag- 
nitude of TpsHP position-dependent. It is however pos- 
sible to impose an average normalization condition. For 
simplicity, we will choose bias l\~g /2 where j3 — Q and 
a — 1. We define the unit vectors: 



(51) 





ei = sin(6'i) cos((y3i) X 






-Hsin(6'i)sin(v5i)y 


+ cos(0i)z. 




62 = sin (^2) cos ((/?2) X 






+ sin(6l2)sin(¥32)y 


4- cos (02 )z, 




q = cos (i^q) X -1- sin ((pq) 


y- 


The 


energy of the SHP is then 
EsHP ,(e'\ 






+ Y^a(g) [sin2(6li) + sin' 


' {02)] 



(52) 



+ 77^ (9) s"^ (6'i)cos(29?i 
lo 



SV'q) 



+ T^^ (9) ^"^^ (^2) cos {2ip2 - Syjq) 



16 
1 



+— c(g)[cos2(0i)+cos2(02)] 



16 

1 

16 
1 

16 



d (q) sin {9i) cos (^2) cos {(fi — fq) 
d (q) sin (6*2) cos (6*1) cos {ip2 - ^Pq) 



In this equation, 6*1 , tpi and O2 , 'P2 must be choosen to 
make ei • 62 = 0. For values of q where all interactions 
a{q) , ...,d{q) are negative, inspection of Eq. ([52]) and 
Fig. [3] shows that the lowest-energy solution is obtained 
when 9i = 0, 6*2 — 7r/2 and (^2 = 'Pq- The vector 62 is 
free to take any orientation in the x~y plane. The plane 
of polarization of the helix is thus z — e2 and its energy 
is given by 



E 



SHP 

N 



/(- 



11 
32 



1 



:/3A, 



(53) 



B. Helical phase 

The exact helical phase (HP) given by the solution of 
the Hartree-Fock equations is too complex to study an- 
alytically and so we give the numerical solution below. 
However it is instructive, in order to understand the in- 
stability of the UP, to look at the energy of a simple 
helical phase (SHP) with only one wave vector: 



TPSHP (r) 



2S 



po [ei cos (q • r) -I- 62 sin (q • r) 



a 
2S 



rjz, 
(50) 



+ Y^iaiq) + b{q)+c{q)+d{q)). 

The sum of the four interactions is plotted in Fig. |31 
From this figure, we see that the wave vector that mini- 
mizes the energy of the helix is qo£ « 2.3. Note that an 
helix in the x — y plane (with 9i = 02 = 7r/2, ipi = 0,if2 = 
7r/2) has a higher energy given by 



E', 



SHP 

N 



-) 

nil 



11 
32 



|-^/3AB + ^a(g). (54) 



Fig. [3]shows that the Heisenberg exchange and dipolar 
parts of the energy of the helix (the curve labelled a + 
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b + c) has its minimum at g = 0. We conclude that these 
interactions do not favor the formation of a non- uniform 
state. In our system, the DM interaction is responsible 
for the formation of the helix with a finite wavevector. 
We remark that our helix has its wave vector q parallel 
to 62 instead of n = ei x e2 as is often the case. 



2. Peudospin pattern in real space 

Fig. [6] shows the orbital pseudospin field in real space 
at different biases. The following symmetry relations 
hold 

(Pxix))^^ = {Pxi-x))^(i)_^^, (55) 

(p^(-^))ab = "(^^(~2;))^(i)_^ • (56) 



1. Energy 

When we solve for the exact helical phase using Eq. 
(|15p . we find a multicomponcnt helix with a finite value 
of rjz- We choose 62 = x so that the helix wave vector 
q — nqoSc with qq = 27r/A and the pseudospin vector ( p ) 
rotates clockwise in the z—x plane. We find that the opti- 
mal period of the helix is A/£ == 2.74+0.01(Ab-Ab '72)^ 
close to the value we found above for the SHP, and varies 
only slightly with bias. Since £ (A) = 256/^B(T), the 
period of the helix is of the order 220 A for B == 10 T. 

The energy of the numerical solution is independent 
of the direction of the wave vector q in the x — y plane, 
as expected. As for the SHP, the plane of polarization is 
given by (ei, 82) = (q,z) .The energy of the HP is plotted 
in Fig. [5] along with that of the SHP. The energy of the 
HP is obviously very close to that of the SHP. We find 
indeed that the number of Fourier components in the HP 
is small. 
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In the numerical solution, {p^ (q = 0)) — {py (q = 0)) — 

but {p2 (q = 0)) 7^ (except at the special bias A^ /2). 
As expected, the helix has a finite value for the term rj 
in Eq. ([50)1 because the bias tilts the pseudospin vector 
away from the z axis. 



0.6 




FIG. 6: (Color online) Density profile for the pseudospin fields 
{px (X)) and {pz (X)) in the helical phase for two different 
biases. 

In the hehcal phase (po,o (q)) + (pi,i (q)) — 1 and so 
the orbital density (p„ (r)) = 1/2tt£'^ is constant in real 
space. The densities 



(Po,o (r)) 



(p. (r)) 



(PmW) = 2-(p.(r)). 



(57) 
(58) 
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FIG. 5: (Color online) Energy of uniform (UP), exact helix 
(HP), and simple helix (SHP) as a function of bias at _B = 10 
T. 



Band structure 



For a modulation of the pseudospin density in only 
one direction such as in the helical state, Eq. (ITSI) for 
the single-particle Green's function can be written more 
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simply as 

k 

(59) 
with n, TO, fc = 0, 1 where 

Gn,m {X, t) = - (TrCn^x (t) C^,, (0)) , (60) 

and 

A„,™(X) = -^Ab(5„,„(5„,i (61) 

The order parameters are given by 

{Pn,m (q)> = -^^,„,0 E ^"''^'' (^".™ (^)) ' (62) 



A^c 



X 



with 



{Pn,m (X)) = Grn,n {X , T = 0") . (63) 

We can solve formahy for the Green's function to find 

^n.m (X ) 



Gn,m (X, iuijj 



iojn + p/fi- E+{X)/h 

Bn,m (X) 

iujn + fi/n-E^ {x)/n' 



(64) 



The band structure consists of two bands with dispersion 
E± {X) given by 



E±{X) - i[Ao,o(^) + Aia(X)] 



(65) 



with 



± V[Ao,o {X) - Ai,i {Xf + 4 |Aoa {X)\' 



^n,™W = ^^^tItS^^Wt^. (66) 



Bn,m (X) 



E+ {X) - E^ (X) 

Kjn.n {X) - E^ (X) Sn, 



E+ (X) - E^ {X) ' 
and rl = 1 - n , etc. At T == K, 

(P„,™ {X)) = Brn,n {X) . (67) 

It is easy to show from Eq. (|59l) that, at T = K, we 
have the sum rule 



\{-fix))f^^. 



(68) 



Moreover, From Eqs. (|66l67p . we have also 

(po,o(^)) + (pi,i(^))=l, (69) 



so that the total guiding-center density is unmodulated 
in the spiral phase. Note that the real density n (r) given 
by Eq. ([25|) is modulated however. The modulus of the 
pseudospin vector \{p {X))\ = 1/2 is constant in space. 
This is not the case when there is a two-dimensional mod- 
ulation of the pseudospin texture as we will see in the 
crystal phase. 

The band structure of the HP is shown in Fig. [7] for 
different values of the bias. The lowest band E- (X) is 
completely filled so that the system is insulating. There 
is a continuum of electron-hole excitations in the energy 
range Eeh & [0.56,0.80] {e^ / nt) which is roughly inde- 
pendent of the bias. The band structure has the symme- 
try 



i?+(X)U^=2AW_£;_(_X)|^a 



-Af 



(70) 



where A^-^ was defined previously in Eq. pS)) . 
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FIG. 7; (Color online) Band structure E± [X) in the helical 
state for different biases and for the uniform phase (straight 
lines) at As = or Ab = A^^'. 



4-. Density of states 

We compute the density of states (DOS) from the re- 
tarded single-particle Green's function i.e. 

g{E,^B) = -^Y.^ [Gl^ (q = 0, E)] , (71) 

n 

- -^EE^[Gn.«(^'^'^)]- 
n X 

The DOS is represented in Fig. |8]for A^ = A^^V2- At 
this bias, the lower energy band is the mirror image (with 



respect to a mirror line at the energy A^') of the high 
energy band. More generahy, because of the symmetry 
of the band stucture, the DOS has the corresponding 
symmetry 



9{E)\^^ = g{2^%^ -E 



A^J'-A, 



(72) 



The extrema in the band structure shown in Fig. [7] lead 
to distinctive van-Hove singularities in the DOS as seen 
in Fig. H 
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FIG. 8: Density of states in the helical state at bias As = 
A«/2. 
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in the Brillouin zone. We remark that, in order to capture 
the electron-hole continuum from x^ a computed in the 
GRPA, we must sum over all the reciprocal lattice vectors 
i.e. compute 



Xa,a (k, i^) 



G 



X^,(k + G,k+G,c.) 



(73) 



This function is shown in Fig. [TUl where the electron-hole 
continuum together with the collective modes are clearly 
visible. Note that all response functions are coupled in 
the GRPA equations. Consequently, they all share the 
same poles. However, the weight of a given pole depends 
on the nature of the underlying mode and is not the 
same in all response functions. Electron-hole excitations 
appear as very localized excitations and are captured in 
the response functions at finite G. 
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5. Response functions and collective modes 

We compute numerically the 16 retarded re- 
sponse functions in the GRPA using Eq. (1511) . 
In the non-uniform phases, one must calculate 
Xi',i,j,j (k -I- G, k -I- G, w) for all recriprocal lattice vector 
G. Using the transformations given by Eq. (P^ . we ob- 
tain the response function for the pseudospin operators 

i-e- Xa,b^^^^O,^f' = Pn,Px,Py,Pz- 

We show in Fig. IH] the imaginary part of the re- 
sponse function x^',!^p„ (k, k, w) defined in Eq. ([32]) . This 
function corresponds to the single-bubble approximation 
and does not capture the collective modes but only the 
particle-hole continuum. The continuum appears in the 
range Eeh G [0.56, 0.80] in accordance with the band 
structure calculation. 

The collective modes can be obtained from the poles of 
the imaginary part of the full GRPA response functions 
X^Q (k, k, a;) . To get the dispersion relation, we follow 
the frequencies of these poles as the wave vector is varied 



FIG. 9: Imaginary part of the density response function (o = 
p„) evaluated in the Hartree-Fock approximation (single- 
bubble approximation) in the helical phase for k — (0, 1) 2tt/X 
and As/ (eV^^) = 1-0. 

Fig. [11] shows the dispersion of the first collective 
modes of the helical phase for As — 1.0 along the di- 
rection of the pseudospin modulation i.e. along kx ■ (The 
absence of points in the dispersion of some of the modes 
is a numerical artefact.) The real-space pattern is pe- 
riodic with period A along k^ but there is no period- 
icity in the dispersion in the ky direction as shown in 
Fig. [12] We have indicated in Fig. [12] the region of 
the electron-hole continuum where the collective modes 
are damped. The higher-energy collective modes are less 
dispersive and correspond to more localized excitations 
which eventually vanish in the electron-hole continuum. 
Since we are dealing with a continuous structure, the 
number of collective modes is not finite. This is also true 
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FIG. 10: Imaginary part of the density response function 
(a — pn) evaluated in the GRPA in the hehcal phase for 
k = (0, 1) 27r/A and As/ (eV W) = 1-0. 



in the Skyrmion crystal phase that we study in the next 
section. 
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FIG. 11: Dispersion relation of the first collective modes of 
the helical phase along k^ for As/ie^ /ni) = 1.0. The rotation 
of the pseudospins is along the x direction with a period A. 



The unusual oscillations in the dispersion relation in 
the ky direction are due to the phase factor 7q.q' = 
exp (-iq X q'^V2) = exp {-i (k + G) x (k + G') >V2) 



in Eq. 



In the spiral phase, G = (27rn/A) x and 



FIG. 12: Dispersion relation of the collective modes of the 
helical phase along ky for Ab / {e^ / k,£) = 1.0. The rotation of 
the pseudospins are along the x direction with a period A. The 
two lines indicates the region of the electron-hole continuum 
where the collective modes are damped. 



k is a vector in the first Brillouin zone. It follows 
that 7q,q/ = exp {—iky {Gx — G'^.) ^^/2) . In the direction 
ky — 0, the phase factor 7q,q' is constant and there is 
no oscillation. Along /cj,, however, the phase factor is 
periodic with periods X^ /mnl'^ in units of 27r/A. In Fig. 
fT2l we can clearly see the periods with m — 1,2,4 using 
A/£«2.74. 

The gapless mode of the helical phase is a phonon mode 
corresponding to a global translation of the density n (r) 
in a space (for q = 0) accompanied by an in-phase z — x 
rotation of the all the pseudospins. As can be seen from 
Eq. ([55)1 , this is a global symmetry of the Hamiltonian. 
At small wave vector k, the dispersion of the phonon 
mode has uj{kx,ky =0) ~ k^ (the wave vector of the 
helix is q =qo^) while cu {k^ = 0, ky) ~ ky-^. 

The gapless x — y pseudospin mode of the uniform 
phase acquires a gap in the helical state. Since the bias 
fixes the z component of the pseudospin in the uniform 
phase, this mode corresponds to an oscillation of the 
pseudospins at pz (r) constant when k = 0. In the he- 
lical phase, the mode that corresponds to this motion is 
gapped because the DM term sets a preferred plane of 
rotation for the pseudospins i.e. the direction of n and 
q in Eq. ([5T|) are related. We remark that in a dou- 
ble quantum well systems where a stripe state occurs at 
V = I'm. higher Landau levels, both the phonon and pseu- 
dospin modes are gapless in the absence of tunneling—. 
There is no DM term in the Hamiltonian of the stripe 
state and both motions correspond to a symmetry of 
the Hamiltonian. A plot (not shown) of the susceptibili- 
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ties Xaa C* ~^ 0,uj) shows that the lowest-energy gapped 
mode in Fig. [T2]has a weight in only xf^,p^ and xf ,p ■ It 
thus seems likely that this mode is the gapped pseudospin 
X — y mode. 



C. Skyrmion crystal phase 

The skyrmion crystal phase occurs on both sides of 
the helical state in the phase diagram. Fig. [T3l shows the 
pseudospin fields defined in the guiding-center represen- 
tation by Eq. ^ at bias As = 0.2. From Eq. dM]), the 
parallel component of ( p (r)) is directly related to the 



(1) 



physical electric dipoles. The crystal at A^ 
electron-hole conjugate of that at A^ i.e. 



n(r)| 



A>„^'-Af 



(P^(r)>L(i) 



A^„'-Ab 



2-«(r)lA«, 
-(p.(r))lA«- 



As is the 

(74) 
(75) 



The pseudospin vorticity in the x — y plane is however 
the same for both crystals. As in the helical state, the 
electron-hole conjugation applies to the other properties 
described in this section. 

In the crystal phase, the pseudospin ( p (r)) is not 
constant in space. In fact, ( p (r)) is a pseudospin den- 
sity not a unit field. In Fig. [T31 the pseudospin den- 
sity has not been normalized and so {pz (r)) is greater 
than 1/2 in some regions of the crystal. The guiding- 
center density {p (r)) = 1 is however constant. The 
density in each orbital is given by Eqs. (|57I58I) . The 
difference between the orbital and crystal cases is that 
{pz (r)) e [—1/2, 1/2] in the orbital state while {pz (r)) is 
positive for A^ < A^Y^ and negative for As > A^72 
in the crystal phase, in the regions where the crystal is 
the ground state. 

The crystal state is constructed by assuming that the 
number of vortices is equal to the number of electrons. 
Since the number of flux quanta is equal to the number 
of electrons at filling factor j/ = 1, we have v — 2Trn£'^ 
with 71 = 1/ea^ where e = a/3/2 for a triangular lattice 
and the lattice constant is given by a/i — \/2t^ je . The 
lattice spacing is thus controlled by the magnetic field. 

The pseudospin pattern at each crystal site in Fig. [T3] 
ressembles that of a meron (half a Skyrmion) with neg- 
ative charge. The z component of the pseudospin is up 
at the center and goes in the x — y plane in-between the 
crystal sites. The vorticity of the field ("^|| (r)) is posi- 
tive. Evaluation of the real density n (r) shows that each 
meron is placed on a uniform background of negative 
charge n (r) i.e. the density does not go to zero between 
two crystal sites. Approximatively half of the real density 
is in the uniform background and half is in the merons. 
We could thus describe the state at As < A^y2 as 
a crystal of charged merons on top of a uniform back- 
ground of electrons. For A^ > A^Y2, the vorticity 
is still positive, but the z component of the pseudospin 



is down at the center of the vortices and we have posi- 
tively charged anti-merons (or holes in the electronic den- 
sity). Note that the electronic density inside each meron 
is not quantized in our approach since we work with a 
pseudospin density and the pseudospin modulus changes 
with position. We could think of the crystal state has a 
two-dimensional charge-density-wave with an amplitude 
that can change continuously with bias. We use the term 
"Skyrmion crystal" in a loose way to refer to that state. 




FIG. 13: (Color online) Pseudospin pattern l^px (r)) , {py (r)) 
in the Skyrmion crystal at bias As/ (e^ /ntj = 0.2. 



1. Density of states 

The density of state is shown in Fig. [M] It is closer 
in shape to the DOS of the helical phase than to that of 
a crystal as can be seen by comparing Fig. [T3] with the 
DOS of a Skyrmion crystal dX v — 1.2 and Ag = 1.28 
shown in Fig. 13 of Ref. il2i . 



2. Collective modes 

The dispersion relation of the collective modes of the 
skyrmion crystal is shown in Fig. 1151 Only the first 
low-energy modes are shown. The modes become more 
dense at higher energy until the electron-hole continuum 
is reached. From Fig. I14[ the continuum is in the range 
Eeh e [0.47, 0.90] . The dispersion for A^^^ - 1\b (not 
shown) is exactly the same as that for A^ as expected. 

Since our calculation does not include disorder, the 
lowest-energy mode is a gapless phonon mode as in the 
helical state. For As = 0.1, the dispersion is isotropic 
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FIG. 14: Density of states of the Skyrmion crystal phase at 
bias Ab/ (cVk/) = 0.2. 



is the gapped pseudospin x — y mode. 

The energy difference between the crystal and helical 
phases is quite small, of the order of a few percent. It 
may be then, that these phases will be very sensitive to 
disorder as well as thermal and quantum fluctuations. An 
evaluation of the effect of these perturbations is however 
beyond the scope of the present paper. 



VI. ELECTROMAGNETIC ABSORPTION 

The collective modes of the helical and Skyrmion crys- 
tal phases can be detected in electromagnetic absorption 
experiments. With k = h for graphene on Si02 sub- 
strate, we have v^ = e^ jl-nhKi = 0.43\/^ THz with B 
in Tesla. The frequency of the collective modes at q = 
in both the helical and crystal phases are in the range 
V e [0.1, 0.6] z/Q « [0.14, 0.84] THz for B = 10 T. 

Theoretically, the absorption can be related to the 
current-current correlation function as explained in Ref. 
[g. We give here only the main results of this calculation. 

We write the current operator, projected onto iV = 
and valley K' as 



dH%^, 



dA'' 



(76) 




0.5 1 
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FIG. 15: Dispersion relation of the first low-energy collective 
modes of the Skyrmion crystal at As/ {e^ / K,tj = 0.2 along 
the path F — J — X — Fin the irreducible Brillouin zone. 



with w ~ k}'^ at small wave vector as in a Wigner crystal. 
The pseudospin x — y mode which is usually gapless in a 
Skyrme crystaUi is now gapped because of the DM term 
in the pseudospin Hamiltonian. A plot (not shown) of 
the susceptibilities x^a C^ ~^ 0, lu) shows that the lowest- 
energy gapped mode in Fig. [15] has substantial weights 
in X 



p^,p^ 



and X 



Py^Py 



It thus seems likely that this mode 



where A*^ is the vector potential of the external electro- 
magnetic field and H^, was defined in Sec. II. We de- 
fine the current-current correlation function Matsubara 
Green's functions as 



XJ.,jA^) = -j{TrJc.{T)J0m 



(77) 



71 



47rm*2£4Xp^,P^(^)> 



where A is the area of the 2DEG and 

Xp^.P, (t) = - {Trp^ (0, r) pp (0, 0)) (78) 

with a, f3 — X, y. Note that a, /3 are defined hyx = y,y = 

X. 

The electromagnetic absorption for an electric field ori- 
ented along the direction a is given by 



Pc.{i0) 



ic 



xl,j„ (w) 



E'. 



(79) 



-KxU^V^^o^ 



xy 



71 



X%,p^ (^) 
to + i5 



where we have assumed a uniform electric field E {t) = 
Eoae~^'^* with polarization a and taken the analytic 
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continuation if2„ — > uj + i5 oi XJa,Ja (Oi*^n) to get 
the retarded response function. The response functions 
^?q pb ('^) ^^^ calculated in units of h/ (c^/k^) so that 
Pa (w) is the power absorbed per unit area. In Eq. ([7^ 
we have neglected a diamagnetic contribution to the cur- 
rent response which becomes important at low frequen- 
cies. 

The absorption due to all but the gapless mode in the 
helical phase is shown in Fig. [121 Comparing with Fig. 
[T2l we see that all the other modes with the exception 
of the third mode (with frequency near 0.3 [e^ /HkVj) are 
optically active. The lowest gapped mode is the most 
intense one and its excitation is strongly sensitive to the 
orientation of the polarization vector of the electromag- 
netic wave. This is true at all bias voltages. The absorp- 
tion frequency does not change significantly with bias. 
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FIG. 16: Electromagnetic absorption Pa (a;) in the helical 
phase for polarization in the x and y directions at bias 
Ab/ {e^/Kt, = 1.0. 

The absorption spectrum for the Skyrmion crystal 
phase is shown in Fig. [171 In this case, all the modes (ex- 
cept the gapless phonon mode) are equally active in the 
absorption and the absorption does not seem sensitive to 
the polarization. There are thus qualitative differences 
between the absorption in the helical and crystal phases 
that should help to observe the transition between these 
two phases. Note that in the UP the orbital pseudospin 
mode is gapless and does not lead to absorption. 



VII. DISCUSSION 

We have studied the phase diagram of the 2DEG in a 
graphene bilayer in the Bernal stacking at filling factors 
1/ = —1, 3 in Landau level A^ = when orbital coherence 
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FIG. 17: Electromagnetic absorption Pa (a;) in the crys- 
tal phase for polarisation in the x and y directions at bias 
As/ {e^lntj = 0.2. 



is present in one of the layer. Our model uses a tight- 
binding Hamiltonian that is simplified by working in the 
two-band model. This simplification is justified since we 
are interested in the low-energy excitations of the 2DEG. 
Moreover, we neglected the warping term, an approxima- 
tion that is valid at sufHciently strong magnetic field. Fi- 
nally, we assumed complete spin polarization. This last 
approximation may fail i^ = — 1 at finite bias when levels 
with spin up mix with levels with spin down. A more 
exhaustive studyS shows that, when this approximation 
is not made, it is still possible to find states with orbital 
polarization although at different filling factors and for 
different ranges of bias than those studied in this paper. 

The physics of the orbital coherent state is due to a 
competition between the Coulomb exchange interaction 
and a Dzyaloshinskii-Moriya interaction between the or- 
bital pseudospins. This competition is responsible for a 
phase diagram where the ground state evolves from a uni- 
form state with collectively oriented orbital pseudospins 
at small bias into a Skyrmion crystal state and then into 
an helical state where the pseudospins rotate in space. 
If the bias is further increased, the helical states transits 
into the Skyrmion state again and then back to the uni- 
form state. All three states can be distinguished from 
their density of states and collective excitations. 

As was shown previously^, the Goldstonc mode due to 
spontaneous orbital coherence in the uniform phase has 
the peculiarity of being highly anisotropic. The disper- 
sion is still highly anisotropic in the helical state which 
is modulated in one direction only but is isotropic in the 
Skyrmion crystal state. The three phases have one gap- 
less mode and several gapped modes. We have calcu- 
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lated that these latter modes lead to absorption in the 
far-infrared region of the electromagnetic spectrum. In 
the helical state, the absorption intensity is very sensitive 
to the orientation of the polarization vector of the elec- 
tromagnetic wave. This is not the case in the Skyrmion 
crystal phase. 

The helical and Skyrmion crystal phases each support 
a gapless phonon mode which is accompanied by a mo- 
tion of the pseudospins. We think that, in the presence 
of disorder, these Goldstone modes should lead to strong 
absorption of electromagnetic waves at very small fre- 
quencies. This calculation is beyond the scope of this 



paper and we leave it for further work. 
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